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ABSTRACT 

Eighty-two ocean bottom gravity stations in southern 
Monterey Bay were occupied in the summer and fall of 1972 from 
the R/V ACANIA. A land gravity survey of ten stations about 
the perimeter of the Bay was conducted in the spring of 1972. 
Gravimeters employed were LaCoste and Romberg Models H6G and 
G-17B, respectively. 

Conventional steps in data reduction are discussed, and a 
terrain correction theory unique to ocean bottom gravimetry 
is presented. The complete Bouguer anomaly (CBA) field for 
bottom and shoreline surveys is included. 

The geological interpretation of the gravity data is 
discussed briefly. Sub-bottom structure of southern Monterey 
Bay as determined by seismic reflection is verified by the 
CBA field, and a calculated density contrast between the 
basement granodiorite and overlying sedimentary strata is 
found to be realistic. The data supports the existence of a 
fault oriented beneath the Monterey Submarine Canyon. 
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INTRODUCTION 



I . 

V 

Geological oceanographers have three principal means of 
investigating structures below the ocean floor: seismic, 

magnetic, and gravimetric surveys, the latter of which is the 
subject of this writing. The relatively small variations in 
the acceleration of gravity over the surface of the earth are 
largely due to such obvious factors as position, elevation, 
crustal density, and local terrain; once these effects which 
may be unrelated to substructure irregularities are removed, 
however, the investigator is left with a meaningful residual 
in the form of the complete Bouguer anomaly (CBA) . This mea- 
sure of anomalous attractive force per unit mass then consti- 
tutes a basis for the interpretation of the upper few hundreds 
or thousands of feet of the earth's heterogeneous crust. 

This study is based upon measurements made at 10 land 
stations and 82 underwater stations which were occupied in the 
spring and summer of 1972. The bottom gravity data collected 
was obtained aboard the Naval Postgraduate School's (NPS) re- 
search vessel R/V ACANIA (Fig. 1) . 

A. SURVEY PURPOSE 

The main purpose for conducting this survey was to learn 
more about the structural geology of southern Monterey Bay by 
obtaining gravity data in areas too shallow for larger vessels 
capable of surface gravimetry. Such vessels must stabilize 
their course and speed for several minutes after each turn 
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Figure 1. U. S. Naval Postgraduate School's Oceanographic Research 
Vessel R/V ACANIA. 




before useful data can be obtained; this usually places pro- 
hibitive restrictions on how close to shore the vessel may 
work, depending upon local bathymetry. But such geographical 
restrictions do not apply to bottom gravimetry, since the 
research vessel can slowly decelerate to even very shallow 
stations and stop for meter lowering and measurement. The 
R/V ACANIA was able to survey within 100 yards of the beach in 
water less than 50 ft deep. Additionally, bottom gravimetry 
is generally more accurate than sea surface gravimetry due to 
the stability of the measuring platform, the sea floor. An- 
other advantage of bottom gravimetry is its closer proximity 
to the density contrasts that produce a gravity anomaly; sur- 
face ship measurements must contend with the attenuating effect 
of vertical distance from the meter to the substructure (Beyer, 
et al . , 1966). This means that the deeper the station, the 
greater the accuracy of bottom gravimetry over surface gravi- 
metry (assuming equal meter precision) . 

After a brief area orientation, this report will describe 
some basic theory of gravity measurements, specific details on 
equipment used, and various procedures employed during the 
survey itself. This will be followed by a section explaining 
each step in the reduction of the raw data, and an analysis of 
the results in terms of geological significance will conclude 
the report. 

B. AREA DESCRIPTION AND GEOLOGICAL SETTING 

Monterey Bay is located about 70 miles south-southeast of 
San Francisco, California. It is characterized by gently 



13 



sloping sandy bottoms in the north and south, divided by one 
of the largest submarine canyons in the world. The area 
encompassed by the author's survey includes that portion of 
southern Monterey Bay within the 50- fm depth contour from 
Pt. Pinos to Moss Landing (Fig. 2) . 

With the exception of the rocky bottom in the extreme 
southwestern portion of the survey area and the predominance 
of green muds near the Monterey Submarine Canyon, a sandy 
ocean floor exists in southern Monterey Bay. Cretaceous 
Santa Lucia granodiori te , readily visible at Pt. Pinos and 
constituting much of the adjacent sea floor, is buried to in- 
creasing depths beneath less dense sediments and sedimentary 
rocks as one proceeds toward the Salinas River mouth-Moss 
Landing area (Fig. 3) . Much of the overlying material to the 
north and east is the Monterey Formation, a Miocene rock 
consisting largely of silicious mudstone, diatomite, and marine 
shale, interbedded with beds of opaline chert. Above the 
Monterey Formation are thin layers of the Pliocene Paso-Robles 
Formation and Pleistocene Aromas red sands. Thinner layers of 
deltaic deposits from the Salinas River have been laid down 
above these formations during the Holocene period. Additionally, 
the Pliocene Purisima Formation, abundant in northern Monterey 
Bay, has recently been found along the south wall of the 
Monterey Submarine Canyon (H. G. Greene, personal communication, 
1973) . 
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Figure 2. Survey Area. 




Figure 3. Geological Map of the Survey Area 
(after Greene, 19 70). 



16 



C. PREVIOUS WORK 



Past geological investigations of Monterey Bay have 
mainly concerned themselves with the origin and present 
structure of the Monterey Submarine Canyon, as seen in the 
studies by Shepard (1948) , Martin (1964) , and Martin and 
Emery (1967) . Geological research in the shallow water 
regions of the Bay has been conducted, but most of these 
efforts investigated specific geological properties, such as 
bottom sediment type and distribution. The first extensive 
and detailed report of the geology of southern Monterey Bay 
was recently presented by Greene (1970). This work reviewed 
the results of previous bathymetric data, rock dredges, and 
grab and core samples, as well as the author's own seismic 
reflection survey results, in order to make inferences con- 
cerning the salt water intrusion problem in the aquifers of 
the lower Salinas Valley. 

While gravity data has heretofore been totally absent in 
the shallow off-shore regions , various investigators have 
surveyed adjacent land areas. Fairborn (1963) , Sieck (1964) , 
and Ivey (1969) have surveyed the local Monterey-Salinas -Ft . 
Ord region, and Bishop and Chapman (1967) compiled land 
gravity data from many sources in preparing their Bouguer 
gravity map of the area. 
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II. THEORY OF GRAVITY MEASUREMENTS 



A. INSTRUMENTS 

When the variation of a quantity is extremely small com- 
pared to the value of the quantity itself, measurement of that 
variation requires some method of exaggeration in order to 
achieve desired sensitivity. Such is the case at hand: 
gravity values at the equator and the poles are approximately 
978 and 983 gal, respectively (owing mainly to the poleward 
decreases of the local radius of the planet and the centrif- 
ugal force due to its rotation) ; hence, the entire range of 
possible variation relative to the mean value is about 0.5%. 
Additionally, this small total range of values is about a 
hundred times greater than the range observed in Monterey Bay, 
so the requirement for precise instrumentation is clear. Most 
of the numerous different gravimeter designs that have become 
operational over the years illustrate this idea. Depending 
upon operating principle, they fall under one of two general 
categories: stable and unstable systems. 

1 . Stable Gravimeters 

Stable gravimeters use a spring-balance system wherein 
the attractive force on a known mass suspended from the end of 
a spring is linearly proportional to the resultant elongation. 
In the case of the Hartley gravimeter, the minute elongation 
is made more measurable by passing a small, horizontally- 
hinged beam through the spring-mass system and attaching a 
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mirror close to the unhinged end, where beam displacement is 
greatest (Fig. 4) . Light from a simple remote source is re- 
flected onto a graduated scale for the observation of relative 
differences. A micrometer screw attached to an auxiliary 
spring is located near the hinged end of the beam so that the 
light beam may be brought back to some fixed reference position 
The amount the micrometer screw is turned to accomplish this 
zeroing process is thus a measure of the difference between 
the local gravity value and the value associated with the re- 
ference position. The system is termed stable because the only 
force acting to counter gravity is the spring restoring force. 

2 . Unstable Gravimeters 

For deflections of the mass in the unstable (or astatic 
balance) gravimeter, the force of gravity is supplemented in 
increasing the spring elongation's restoring force by a third 
force which acts in the same sense as that of gravity itself. 
The Thyssen gravimeter is a typical example (Fig. 5) . The 
clockwise moment of an auxiliary mass above the pivot point 
results in additional beam displacement which can be kept 
linearly proportional to gravity. Thus, greater sensitivity 
is the chief advantage of the unstable gravimeter; for a given 
change in gravity, the unstable device will undergo greater 
beam displacement than that of a comparable stable gravimeter. 

B. OCEANIC VERSUS LAND GRAVIMETRY 

In the case of land gravity meters, systems like those 
described above are enclosed in an instrument roughly the size 
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Figure 4. Schematic of a Stable Gravimeter: 

the Hartley Model (after Dobrin, 1960). 




Figure 5. Schematic of an Unstable Gravimeter: 

the Thyssen Model (after Dobrin, 1960). 
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of a kitchen toaster, and are therefore as mobile as the 
operator. The mass is manually clamped in place after each 
reading and the only electrical requirement is a battery to 
provide power for the light source, scale illumination, and 
temperature regulation (necessitated by the temperature- 
dependence of the expansion coefficients of the springs) . 

With regard to oceanic gravimetry, however, the situation 
is vastly more complicated. Three platform options are: 
surface ship, submarine, and bottom measurements; each has 
its own drawback. Surface ship gravimetry requires intricate 
sensors to measure the effects of ship-motion accelerations. 
Submarines designed for oceanographic research are still quite 
rare; military submarines have severe space limitations and 
the more pressing function of national defense to attend to. 
Bottom gravimetry (conducted with operators accompanying the 
meter to the bottom in diving bells prior to the 1950 's) in- 
volves the expensive business of remote control. In addition 
to these monetary problems, all types of oceanic gravimetry 
must contend with the formidable problem of accurate navigation, 
a difficulty not common to most land surveys. 

Other difficulties are inherent and peculiar to bottom 
gravimetry. The gravimeter itself must be enclosed within 
some type of waterproof chamber strong enough to withstand 
pressures of many atmospheres. Also, a large number of 
separate electrical conductors, normally within the insulation 
and armor of an oceanographic cable, must be employed to effect 
the numerous recordings and operating modes. This requires a 
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specialized winch with a multitude of slip-rings. Rocky or 
steep bottoms can incline the instrument to an extent beyond 
the range of the leveling gimbals. But perhaps the most 
challenging problem of all is the requirement that the re- 
search vessel maintain its position above the stationary 
meter, despite drift induced by wind and current. If this 
requirement is not met, and the winch operator has not lowered 
sufficient excess cable, the meter may be dragged and jarred 
while in the unclamped mode, thereby damaging the meter's 
delicate spring suspension system. 
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III. EQUIPMENT 



A. GRAVIMETERS 

Of the two general types of meters discussed in the pre- 
vious section, the LaCoste and Romberg models fall in the un- 
stable category. Figure 6 is a drawing of the main components. 
When the mass is released from its protective clamps for the 
purpose of taking a reading, its motion to an equilibrium 
position is impeded by very heavy air damping; it therefore 
takes a finite time interval for the mass, support beam, and 
light beam to come to rest. The speed of the motion of the 
light beam is a function of the position of the measuring screw 
and is indicated by a beam position galvanometer which may be 
connected to a strip chart recorder. The slope of this record 
of time versus speed is a measure of how close the measuring 
screw is to the equilibrium position, which, in turn, is a 
measure of gravity. 

The zero-length spring in Figure 6 was first introduced by 
LaCoste in 1934 for use in an innovative long period vertical 
seismograph. Its unique method of manufacture renders elon- 
gations or compressions proportional to increases or decreases 
in gravity itself, resulting in increased sensitivity. 

The calibration of most other types of gravimeters depends 
on weak auxiliary springs which are connected directly to the 
measuring screws (Figs. 4 and 5) . This is avoided in LaCoste 
and Romberg gravimeters by employment of a lever system 
connected to the measuring screw. The light beam is thereby 
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Measuring 




Figure 6. Simplified Diagram of the LaCoste 

and Romberg Gravimeter (after LaCoste, 
1967 ) . 
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brought to the equilibrium position by actually changing the 
location of the zero-length spring's upper support. 

1 . LaCoste and Romberg Model G-17B Geodetic Land Gravimeter 
The LaCoste and Romberg model G-17B geodetic land 

gravimeter used for the shoreline survey is an instrument 
housing the system just described in a 5-3/4 inch x 6-1/2 
inch x 8-1/2 inch case weighing about 5 lb. This model was 
first introduced in 1956, and operates in the field on battery 
power over a 7000 mgal range. Instrument drift is less than 
0.5 mgal/month due to a thermistor-transistor heater control 
system. Sealed against atmospheric pressure changes, the G-17B 
gravimeter is nulled by adjusting the measuring screw while 
observing the light beam positon through a microscope eyepiece. 

2 . LaCoste and Romberg Model H6G Underwater Gravimeter 
The LaCoste and Romberg Model H6G underwater gravimeter 

used for the bottom survey is a complete system which permits 
an instrument like the G-17B to be used on the ocean floor. 
Simply stated, it is a remote-controlled land gravimeter in- 
side a shell of two thick aluminum hemispheres supported by a 
triangular base (Fig. 7 and 8) . A brief description of the 
features of this model is presented in Table I. 

B. AUXILIARY EQUIPMENT 

The bottom gravimeter and all of the associated equipment 
were loaned to the NPS Oceanography Department by the U. S. 

Naval Oceanographic Office in Suitland, Md. Engineered and 
assembled by the Geophysical Equipment Manufacturing Co. of 
Houston, Tex., the auxiliary equipment's operation may be 
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Figure 7. Kodel H6G Gravimeter Ready for Use. 
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Figure 8. Internal View of the yodel 16 G Gravirster 




Table I. Features of the LaCoste and Romberg 
Model H6G Underwater Gravimeter 



Ba se Size 
(with legs) 


39-inch equilateral triangle 


Height 


29 inches 


Weight 


350 lb, including enough weight 
to insure acceptable sink rate 


Range 


World wide (7000 mgal) 


Optimum Accuracy 


+0.02 mgal 


Realistic Accuracy 


+0 . 10 mgal 


Power Source 


115 vac, 60 Hz 


Routine Maximum 
Operating Depth 


600 ft 


Maximum Operating 
Depth Tested 


2900 ft 


Maximum Bottom 
Slope Permissible 


15° 


Depth Indication 


Pressure sensor (0.5% accuracy) 
connected to casing port 


Leveling System 


Optically-controlled servo motors 
driving rack and pinion assemblies 


Flood Indication 


Open circuit, 1/4-inch gap, inside 
bottom of lower hemisphere 
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summarized as follows : a gasoline engine provides mechanical 

energy to a hydraulic pump transforming the rotary motion into 
hydraulic pressure (Fig. 9). This hydraulic pressure is used 
to operate both the winch and A-frame using two, two-way con- 
trol valves. One valve controls a hydraulic motor which con- 
verts the fluid pressure back into rotary motion to drive a 
heavy-duty chain link connected to the winch reel's shaft. 

The A-frame is hinged at deck level and actuated from above by 
two hydraulic pistons controlled by the other valve. 

With the gravimeter connected electrically to the termina- 
tion at one end of the cable, the other end is routed inside 
the hollow shaft of the winch and out the side to a slip-ring 
assembly. The outer side of the slip-ring assembly is wired 
to a cannon plug connection on the control box, which in turn 
is connected through a rectifier and an isolation transformer 
to one of the ship's 115 vac outlets, the basic power source 
for the entire electrical system. 
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Figure 9. Schematic Diagram of the Auxiliary Equipment. 



30 



IV. PROCEDURE 



This survey was part of a larger effort to investigate 
off-shore gravity values over a more extensive area: two co- 

workers (also NPS students) recorded data in northern Monterey 
Bay and Carmel Bay. Since bottom gravity studies are by no 
means a one-man operation, they were assisted by the author 
and provided similar assistance to him during his survey. 

A. GRAVIMETER CALIBRATION 

Prior to beginning the shoreline survey on 26 April 1972, 
the land meter was calibrated by the author and co-workers 
under supervision of U . S. Geological Survey (USGS) personnel. 
Readings were taken at several benchmarks from USGS head- 
quarters in Menlo Park (USGS 1 JD) to Skeggs Point Scenic View 
(U. S. Coast and Geodetic Survey Station B-388) , a frequently 
used calibration route spanning a range of 137.20 mgal (Chapman, 
1966) . Subsequent reduction of counter values to gravity values 
resulted in a difference of 137.15 mgal over the entire range. 
Smaller differences were found at the intermediate benchmarks. 

The underwater gravimeter was initially calibrated by the 
manufacturer prior to shipment to Monterey in the spring of 
1972. After the survey was completed, another calibration 
check of the meter was in order. Early in 1973 the author 
occupied CA-259 (Woollard Airport Base WA-84) at the Monterey 
County Airport and a benchmark located at the base of the steel 
tower on the Monterey Coast Guard pier, designated WH-29 by 
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Woollard and Rose (1963) . The meter readings resulted in a 
gravity difference less than 0.5 mgal greater than the 22.5 
mgal established value (Chapman, 1966) . The difference was 
probably due to the recent construction of a second terminal 
building next to CA-259. Since the airport value is lower 
than the pier value, and since the effect of the new building 
is an additional upward attraction, the difference between the 
two reference stations has probably increased due to the new 
construction. In view of this consideration, as well as the 
fact that the published location of the CA-259 benchmark could 
not be found exactly (due to the construction) , it is estimated 
that the observed gravity difference between the airport and 
pier stations was within 0.1 mgal of the true value. 

B. SHORELINE SURVEY 

A preliminary gravity survey of the perimeter of southern 
Monterey Bay was undertaken in late April, 1972, with the land 
gravimeter discussed above. The information from this shore- 
line survey of 10 stations was later used to verify the near- 
by CBA values and isoline trends of the underwater survey, 
since land gravimetry is an order of magnitude more accurate 
than oceanic gravimetry. (While land gravimeters are read to 
the nearest 0.01 mgal and subsequent corrections' accuracies 
can lead to CBA values of 0.1 mgal accuracy, remote-controlled 
bottom gravimeters, read to the nearest 0.1 mgal, can result 
in a CBA accuracy of at best 1 mgal.) In addition to verifying 
the bottom survey results, these land stations also served to 
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tie in the combined survey with previously surveyed areas in- 
land. All 10 land stations were located (vertically) within 
a few feet of sea level so that the small magnitude of the 
subsequent gravity corrections would yield better CBA 
accuracies. The stations were fairly evenly spaced from Pt . 

Joe to Moss Landing (Fig. 10) . Visual estimates of station 
level relative to the average existing sea surface elevation 
were recorded at each station, so that station elevations rel- 
ative to mean sea level could be calculated later from tidal 
information. Exact locations and depths/elevations of both 
underwater and shoreline stations are included in Appendix A. 

C. SHIPBOARD INSTALLATION AND SEA TRIALS 

When the auxiliary equipment arrived in Monterey in late 
June, 1972, preliminary inspection indicated that it had been 
some time since the system was used. A re-termination of the 
cable and an engine overhaul were required before the survey 
could begin . 

An initial concern was the question of location of the 
equipment aboard the ship. It was decided that the after 
section of the upper level would be most suitable in view of 
the physical requirements associated with meter lowering and 
the intended length of stay of the equipment aboard. The A- 
frame, supporting plate, and gravimeter were stationed on the 
starboard side so that the data recorder in the dry lab directly 
below would be able to observe cable speed and direction of 
motion (Fig. 11). Also, the crane's location in close proximity 
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Figure 10. Station Locations. 
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Figure 11. Auxiliary Equipment Installed Aboard the R/V ACANIA. 



on the starboard side seemed compatible with the plan to bring 
the meter to and from the main deck for each station. The 
winch was placed abeam the A-frame area but on the port side, 
in an attempt to spread out the 4 ton weight of the equipment. 
The engine and hydraulics stand were located on the ship's 
centerline and slightly forward of the other components. 

A brief sea trial period was arranged prior to the start 
of the actual survey in order to acquire a familiarity with 
conditions to be encountered. A trial station was occupied 
west of Seaside, California, in about 100 ft of water. SCUBA 
divers accompanied the meter to the bottom to observe the 
sink rate, take some underwater photographs, and determine if 
excessive cable would foul the feet of the gravimeter. How- 
ever, a strong wind and too little excess cable resulted in 
dragging the meter across the sandy bottom while the meter 
was in the unclamped mode. It was decided to return to the 
ACANIA's mooring to check for possible damage to the mass 
suspension system. Fortunately the counter reading reproduced 
the value obtained just prior to departure for the station, a 
sign that no damage had been sustained. 

This preliminary venture into the Bay showed that the 
stations would have to be occupied with the ship's bow into 
the wind so that the major wind-induced position corrections 
could be accomplished with ahead or astern engine commands. 
Since the ACANIA has no bow-thruster, the bow could be kept 
into the wind with alternate use of one or both of the engines 
and rudder. Also, a feel for the amount of cable payed out as 
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a function of cable angle and tension was obtained, since the 
divers reported that cable fouling about the base of the 
meter was a remote possibility (owing mainly to the cable's 
large diameter and resultant large radius of curvature). 

D. SURVEY OPERATIONS 

Prior to the onset of the survey, selection of station 
locations was governed mainly by the requirement for a fairly 
evenly-spaced grid from which an accurate scalar field of the 
CBA could be extracted. In view of the size of the survey 
area and the expected ACANIA availability periods, it was felt 
that a minimum of 75 stations would be needed to satisfy this 
requirement. A chart showing the spacing of about 100 intended 
stations similar to Figure 10 was provided to the bridge for 
course information. Whether the actual stations occupied coin- 
cided with these intended positions was not of great concern; 
the important thing was to determine the position of each 
station as accurately as possible during and immediately after 
its occupation in order to adjust the planned location of the 
next station for a possible gap in coverage. Stations were 
numbered in order of occupation (except for those on land) . 

Upon completion of the survey, the 92 stations occupied covered 
an area of 50 sq n mile, representing a station density of 
almost 2 per sq n mile (Fig. 10) . 

1 . Navigation 

With the exception of three stations occupied during 
morning fog (and consequently requiring radar positioning). 
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all of the bottom station positions were determined by visual 
lines of bearing from prominent landmarks during hours of both 
daylight and darkness. Three or four such bearings were taken 
at each station, and these were often supplemented by radar 
distances from coastline features. Standard navigational 
procedures were used to estimate the true location of the 
ship's position in the small area where lines and arcs of 
position converged. Most of the three-line position triangles 
were on the order of 0.1 n mile or less on a side, and a good 
portion of the plots were point fixes. Bearing-taking was 
begun as soon as the ACANIA became dead-in-the-water (before 
the meter was lowered to the bottom) and was always completed 
well before the meter was raised after the measurement. When 
it appeared the vessel was drifting, a second set of bearings 
was taken. Differences in meter position and ship's position 
at each station were thus small enough to safely neglect 
consideration. In view of these factors, it is stated with 
a high degree of reliability that station position errors at 
sea were at worst 0.1 n mile. If such an error was in a 
north-south direction the value of theoretical gravity and 
final CBA would differ by only 0.14 mgal . 

Land station navigation was less complex and far more 
accurate. Six of the ten stations were selected for occupation 
because their exact location was clearly shown on the charts 
used. While the other four station locations were transferred 
to the chart from hand drawn maps , it is certain that their 
position accuracy was at least twice as good as that of the 
oceanic stations. 
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2 . 



Measurements 



During the survey the meter was first lowered to the 
sea surface and a depth counter value was obtained by nulling 
the appropriate galvanometer; this was repeated as soon as the 
meter reached the sea floor. Next the leveling mode was 
selected and required flood and tilt checks were made. (Only 
at 2 of the 82 bottom stations did the greater-than-15° bottom 
slope indication appear, requiring repositioning of the meter.) 
After the high-speed leveling switch was activated to make 
coarse leveling adjustments, the meter was put into the read 
mode. In this mode, while fine leveling adjustments continued 
automatically, the beam position and gravity counter switches 
were alternately used to stop the motion of the beam position 
galvanometer needle in the null position. During this process 
the correct counter value was always approached from the same 
direction (namely, from too low a counter value) to avoid 
hysteresis errors . When what was determined to be the correct 
counter value was obtained, an over-adjustment of 0.1 counter 
units was introduced to check for a reversal in the direction 
of beam position motion. This method of obtaining the gravity 
counter value is not only much faster than the slope method 
used with the strip chart record of time versus beam speed, 
but also about as reliable (II. B. Parks, LaCoste & Romberg, 
Inc., personal communication, 1972). After the gravity counter 

value was recorded the mass was again clamped in the deck mode 
and the meter raised. 
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The meter was normally on the bottom for about 4 or 5 
min, 2 or 3 min of which were spent in the read mode. In- 
cluding transit time between stations (but excluding transit 
time to and from the first and last stations of each track) , 
the 82 bottom stations were occupied in a total of 25 hours , 
yielding an average rate of better than 3 stations per hour, 
or about 18 min per station. The deepest stations usually 
required the most time. 

At the shoreline stations the land gravimeter was 
leveled visually with two bubble glasses and manually operated 
screws attached to the frame. The scale was then illuminated 
with a light switch and the mass was manually unclamped. Cor- 
rect counter values were obtained by simultaneously observing 
the beam position through the microscope eyepiece and manually 
adjusting the measuring screw to null it in the reference 
position . 

Fortunately, sny cronization between the counter read- 
ings of the bottom gravimeter and the control box was maintained 
throughout the author's survey. Between the August and Septem- 
ber cruises, however, loss of sny cronization did occur once. 

This was evident when a routine counter reading at Station 1 
was several milligals from its standard value, and the loss 
required opening the inner gravimeter case to read its counter 
value so that the control box counter value could be reset to 
agree . 

An important set of measurements midway through the 
survey, accomplished by repetitive occupation, established a 
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3.20 mgal difference between Station 1 on the floor of 
Monterey Harbor and Station E, the USCG pier benchmark. In 
addition to tieing in all bottom gravity readings with a 
reference benchmark value, this result also related values 
recorded by the two different meters with a station common to 
both surveys. 

3 . Environmental Effects 

From a climatological viewpoint, the summer months are 
best suited, due to wind conditions, to gravity operations in 
southern Monterey Bay. By observing surface bubble relative 
motion as well as cable angle, the ship's master was usually 
able to maneuver the R/V ACANIA to avoid being blown downwind 
of the desired position directly above the meter. The wind- 
induced sea and swell presented more of a problem to survey 
operations by causing noticeable pitch and roll. This occa- 
sionally resulted in the gravimeter hitting the extended A- 
frame during lowering and raising evolutions. 

Readings taken in shallow depths were often charac- 
terized by a slowly fluctuating beam position in addition to 
the normal motion of the beam position prior to obtaining the 
final reading. This required observation of the rate of 
change of the position of the fluctuation's mean so that the 
gravity counter switch could be used to eventually hold this 
mean in a fixed position. Since the automatic averaging 
capability of the control box was inoperative, the averaging 
process had to be done by eye, thereby introducing an esti- 
mated reading error of 0.1 counter units into the final counter 



values . 
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The fluctuating beam position was, in part, the visual 
indication of a form of micro-seismic bottom motion caused by 
surface swell and wave activity, its influence on beam position 
being a function of depth and type of bottom sediment (LaCoste , 
1967) . Additionally, a swaying of the cable was probably felt 
by the meter through the rigidly-attached cable termination. 

The effect of this problem on survey operations was the need 
for more time in the read mode to obtain the correct counter 
value, and, hence, the occupation of fewer stations per hour. 

Fog conditions often prevail in the east central region 
of Monterey Bay during the early morning hours of spring and 
summer. Such was the case during the occupation of Stations 
36, 37, and 38. Its only effect upon the survey was a slight 
deterioration of positioning accuracy; in the absence of 
visual contact with distant landmarks , radar navigation was 
utilized . 
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V. DATA REDUCTION AND PRESENTATION 



The method whereby the raw data collected during the sur- 
vey is transformed into a meaningful CBA picture is discussed 
below. Although much of the theory involved in data reduction 
applies to both land and ocean stations, there are important 
differences. Some of the individual corrections discussed 
apply only to bottom gravity data; a summary of the reduction 
methods applicable to land stations follows in a separate 
section . 

The majority of the required calculations were accomplished 
through use of the NPS IBM-360 computer. Programs were written 
by the author in Fortran language and numerical significance 
was double precision. Copies of these programs are included 
at the end of this report. 

A. OBSERVED GRAVITY 

As pointed out earlier, gravimeters measure gravity differ- 
ences rather than absolute gravity values. Since gravity 
counter units as recorded from the gravimeter control box are 
numerically close to, but not exactly equivalent to, milligals, 
they must first be converted to milligals before addition to 
or subtraction from a reference gravity value. The absolute 
reference used for this survey was 979891.7 mgal , the gravity 
value at WH-29 on the Monterey Coast Guard pier. In practice 
a more convenient intermediate reference station at the mooring 
buoys (Station 1) was used to tie in all of the other station 
readings and to monitor meter drift. 
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The control box counter reading at the pier benchmark was 
3323.05 units; hence, the formula used for computing observed 
gravity at each bottom station was: 

g Q = 979891.7 + [( CV Q - 3323.05 ) x ( 1.03985 )] (1) 

where g Q is in milligals and CV Q represents the observed gravity 
counter value at the station in question. The 1.03985 factor 
converting counter unit differences to milligals w as obtained 
from a table provided by the manufacturer listing (a) a milli- 
gal value associated with each even hundred counter units as 
well as (b) conversion factors for values between the counter 
units listed. Since all of the stations' observed counter 
values fell between 3300 and 3400, only one such conversion 
factor from the listing was required. 

B. LATITUDE CORRECTION 

The first correction (although the order of various cor- 
rections is not defined in any mandatory way) applied to the 
observed value of gravity is normally the subtraction of the 
theoretical value of gravity based on geographical position. 

The determination of theoretical gravity obviously requires a 
precise knowledge of the shape of the earth so that the dis- 
tance from the surface to the earth's center and the centrif- 
ugal force are known for any location. The actual geometry 
has been determined to be very nearly a triaxial spheroid, or 
ellipsoid of revolution having depressions along the two 45° 

latitude lines and a flattening and bulging of the equator 
(Dobrin, 1960). Years of geodetic investigation have yielded 
numerous revised coefficients for the equation of this 
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reference spheriod to more accurately approximate the earth 
geoid. The author's selection of a particular equation was 
governed (as were all aspects of data reduction) by a desire 
to conform with the procedures of USGS and the California 
State Division of Mines and Geology in an effort to render 
the final CBA map compatible with existing land gravity maps. 
The coefficients used were those of the 1930 International 
Spheroid, the theoretical value of gravity being a function of 
latitude only (Dobrin, 1960) : 

g t = 978049.0 (1 + 0.0052884 sin 2 L - 0.0000059 sin 2 2L) (2) 
where L is the latitude. Equation (2) results in a predicted 
south-to-nort.h increase in the value of local theoretical 
gravity of 1.44 mgal/n mile. 

The assumption that the geoid and the reference spheroid 
coincide would appear erroneous at first glance, since, in 
fact, they can differ vertically by as much as 50 m in some 
regions of- the world. However, a more meaningful measure of 
the validity of the effect of this assumption is the change 
in this difference over the area of interest. In the small 
domain of Monterey Bay, the change in the distance from the 
equipotential surface of the geoid to that of the reference 
spheroid is on the order of a few centimeters . It may there- 
fore be concluded that final CBA values will still result in 
an accurate representation of small-to-intermediate scale 
geological features of the upper crust (Grant and West, 1965) . 
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C. EARTH TIDE CORRECTION 



In that the interior of the earth possesses finite 
elasticity, the attractive forces of heavenly bodies (partic- 
ularly the moon and sun) continually act to deform its shape. 

As is the case with the ocean surface, the pertinent manifes- 
tation of this phenomenon is a small scale vertical fluctuation 
of the earth's crust which can affect observed gravity values 
to a measurable extent. The true distance from a point on the 
surface to the earth's center can change by as much as 1 ft 
in a matter of hours, representing a 0.1 mgal change in gravity 
thus, the effect of these earth tides on gravity measurements 
made over such an interval of time must be considered. 

The USGS earth tide Fortran program was utilized to obtain 
the required corrections. Dates and times of station occu- 
pation along wi th a geographically-central location of 36° 

42' N, 121° 52' W were the input parameters. The output values 
were tabulated at 20 min intervals. A copy of this program may 
be obtained from the NPS Department of Oceanography, Code 58Ad. 

D. DRIFT CORRECTION 

Readings taken at Station 1 before and after each track 
were, more often than not, slightly different. These differ- 
ences were attributable to three causes : 

(a) earth tidal variation during the time interval 
between base station occupations; 

(b) ocean tidal variation during the same period; 

(c) meter drift. 
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Since the ACANIA was secured to fixed mooring buoys in an 
east-west orientation, possible north-south variation in the 
actual location of Station 1 before and after each track was 
negligible. To determine true meter drift it is only neces- 
sary to remove the effects of (a) and (b) above from the total 
"apparent drift" actually measured. Recalling that counter 
unit differences are very nearly equivalent to milligal dif- 
ferences, earth tide variation was removed directly from counter 
readings. Similarly, the changes in attraction due to differ- 
ent sea surface tida] heights above the meter were removed. 

(Each 1 ft of water above the meter decreased observed gravity 
by 0.0127 mgal.) The resultant drift values for each track 
were then applied linearly over the range of stations. The 
greatest drift rate was 0.036 mgal/hour. 



E. FREE -AIR CORRECTION 

The formula for the value of the Newtonian portion of 
gravity (in milligals) for a unit mass at sea level (g n ) is: 

G M 



gn = -^2 



(3) 



where G is the universal gravitational constant (6.670 x 
10 - ^cm^g _ ^sec - ^ ) , M is the mass of the earth, and R is the 
distance to the earth's center. The free-air correction 
(FAC) to gravity data accounts for the fact that the measure- 
ment was made at some elevation other than mean sea level, the 
assumed height of the surface of the reference spheroid. (A 
gravity station located exactly at mean sea level, then, would 
not require this correction.) Disregarding the mass between 
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the station and sea level for the time being, the free-air 
correction is obtained by differentiation of (3) (Heiskanen 
and Vening Meinesz, 1958): 

dtJn 

FAC = x (H) = (0.09406 mgal/ft) x (H) (4) 

where II is the elevation of the station above mean sea level. 

A simple modification of this equation yields a formula which 
applies to bottom gravimetry: 

FAC = (0.09406 mgal/ft) x (Water depth - Tide height) 

(5) 

Since equation (3) implies decreasing values with increasing 
distance from the center of the earth, the free-air correction 
is negative for underwater stations. 

A first-order gravity anomaly, the free-air anamoly (FAA) , 
is obtained with the inclusion of this correction along with 
others mentioned thus far: 

FAA = g Q - g t + ET + D - FAC (6) 

where g G , g^ , and FAC are as defined by equations (1) , (2) , 
and (5) , respectively, ET being the earth tide correction, 
and D the drift correction. 

The free-air correction effectively refers station level 
to sea level, and is easily the largest single correction to 
be applied. It is therefore clear from equation (5) that 
accurate values of station depth and tide level are critically 
important if an acceptable degree of precision is to be 
achieved. The pressure sensor in the gravimeter chamber pro- 
vides depth counter indications from the sea surface to the 
sea bottom; the difference in the surface and bottom readings 
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is converted to depth in feet by means of a linear relation 
extracted from sensor calibration graphs provided by LaCoste 
and Romberg. 

All tide heights used were based on the assumption that 
the existing tide level at each station in the survey area 
was identical with that simultaneously recorded in Monterey 
Harbor. Levels were recorded by the NPS tide recording de- 
vice located in the assistant harbormaster's office on 
Municipal Wharf No. 2 in Monterey Harbor. A reference line 
was drawn on these records at 5.5 ft above staff zero. Mean 
sea level is 5.9 ft above Monterey's staff zero (H. V. Maixner, 
personal communication, 1973). Thus, tide height distances 
were measured relative to 0.4 ft above the reference line. 

These tide values were compared with appropriate values from 
the published 1972 tide prediction tables (U. S. Department of 
Commerce, 1972) ; agreement was very good. In order to make 
this comparison it was necessary to determine the vertical 
distance from mean lower low water to mean sea level for 
Monterey, since all tide table values are referenced to mean 
lower low water. This difference was found to be 2.7 ft. 

F. BOUGUER CORRECTION 

While the free-air correction is concerned with changes in 
the denominator of the right-hand side of equation (3) , the 
Bouguer correction (BC) addresses changes in the numerator. 

This correction considers how the attraction between the mass 
within the gravimeter and the earth's center of mass is 
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affected by (a) the presence of material between the station's 
elevation (or depth) and sea level, and (b) the density of 
this material. 

The method of the Bouguer correction first assumes that 
the measurement was made on an infinitely flat ocean bottom, 
so local terrain irregularities are temporarily neglected. 
Removing the upward attraction of an infinite slab of over- 
lying water whose thickness equals the station depth would 
effectively replace that slab of water with air and result in 
a higher value of observed gravity. Since the formula for the 
attraction of such an infinite slab is 2TrGgh, where ^ is 
the density and h is the thickness (Dobrin, 1960) , the first 
part of the Bouguer correction (BCq) , the removal of the sea 
water's attraction, is: 

BC X = 2TTG p w Z (7) 

where Z is the pressure sensor depth and 9 w is the density of 
sea water. The second part of the Bouguer correction (BC 2 ) 
accounts for the increased downward attraction for the sea 
water that would have existed had the measurement been made at 
mean sea level (thus, the previous inclusion of the free-air 
correction is assumed) : 

BC 2 = 2 TT G 9 W (Z-Z t ) ( 8 ) 

where Z t is the tide height relative to mean sea level. The 
third and final part of the Bouguer correction (BC 3 ) considers 
the additional downward attraction that would have existed had 
the material between sea level and the ocean bottom been an 
infinite slab of crustal material instead of water: 
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( 9 ) 



BC 3 = 2Tr G ( P R - 9 W ) (Z-Z t ) 
where p R is the density of average crustal material. Figure 12 
summarizes these considerations schematically. The total 
Bouguer correction is then: 

BC = 2 TP G 9 W Z + 2 Tr G p w (Z-Z t ) + 2 Yt G (?r - p w ) ( Z- Z t ) (10) 

or, combining terms, 

BC = BC X + BC 2 + BC 3 = 2 YT G p w Z + 2Y?G p R (Z-Z t ) . (11) 

The density of sea water was taken as 1.027 g/cm^ , and an aver- 
age crustal density of 2.670 g/cm^ was used for p R . It is clear 
from Figure 12 that the Bouguer correction is positive for under- 
water stations. Once it is included with the other corrections, 
the simple Bouguer anomaly (SBA) emerges: 

SBA = g Q - g t + ET + D - FAC + BC (12) 

Or, from equation (6) , 

SBA = FAA + BC. (13) 

Just as bottom gravity data is reduced to the complete 
Bouguer anomaly in order to compare it with land gravity 
anomalies, its reduction to a gravity anomaly intermediate 
between the free-air and simple Bouguer anomalies will permit 
comparison with FAA values as determined by a surface ship. By 
only adding BC-^ and BC 2 to the free-air anomaly of the bottom 
data, such a value is obtained (Fig. 12-4 applies). This may 
be referred to as a "mass-adjusted" free-air anomaly (FAA 1 ). 

Since the free-air and Bouguer corrections are both func- 
tions of station depth, pressure sensor errors must be exam- 
ined to determine their effect on these corrections. A 1-ft 
depth error produces a combined FAC/BC error of 0.05 mgal 
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(the two corrections are of opposite sign) . Although the 
manufacturer's claim of pressure sensor accuracy is 0.5% 

(this would correspond to + 2 . 3 ft for the deepest station) , 
a more conservative assessment was used. It was assumed that 
the depth determination was in error by 4 ft at most, so the 
combined FAC/BC error is +0.20 mgal or less. 

G. TERRAIN CORRECTION 

Easily the most time-consuming aspect of gravity data re- 
duction is the topographic, or terrain correction (TC) . Often- 
times this correction may be ignored altogether (J. D. 
Rietman, personal communication, 1972) , and this has been 
generally true for previous underwater gravity work (much of 
which seems to have been done in the Gulf of Mexico) . However, 
due to the close proximity to Monterey Bay of very deep abyssal 
plains in the Pacific, as well as such intermediate-scale 
features as the Monterey Submarine Canyon and various coastal 
mountain ranges, this sort of simplification would be unreal- 
istic at best. (In .the final result, quite noticeable differ- 
ences were found in the positions of SBA and CBA isolines 
spaced at an interval of 5 mgal.) Consider a gravity station 
on the ocean floor located on the upper rim of a deep trench. 
The fact that the trench is filled with sea water produces a 
smaller value of observed gravity than would be measured if 
the trench were filled in with solid crustal material. The 
upward attraction of a nearby guyot would also act to decrease 
observed gravity. Thus any terrain deviation above or below 
station level reduces observed gravity relative to its flat- 
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bottom value, both on land and under the sea surface. This 
correction is therefore always positive. 

The generally accepted method for calculating the 
gravitational attraction of topographic irregularities was 
first developed by Hayford and Bowie in 1912. It consists of 
approximating the volume of excess or deficit terrain surround- 
ing the station with a series of concentric cylindrical shells 
of varying height (depending on the local elevation or depth) . 
In plan view this appears as a set of circular zones; the 
zones are lettered alphabetically from the station outward and 
divided into many compartments. This bulls-eve representation 
is then put on glass or acetate templates scaled to conform 
with appropriate topographic or bathymetric charts. By 
centering such a template on the station, the average elevation 
or depth of each compartment is estimated to determine its dis- 
tance above or below station level. This difference between 
compartment and station level is the entering argument for 
tables which have been developed to give the corresponding 
gravitational attraction in milligals. The compartment cor- 
rections (199 in all) are then summed to yield the total 
terrain correction. 

The tables used in this research are modifications by the 
USGS of the work done by Swick in 1942. His tables are based, 
in turn, on Bullard's 1936 modification of the original Hayford 
Bowie paper that presented the combined effects of topography 
and isostatic compensation. These tables assume a removal/fill 
in density of 2.670 g/cm^, and give a 0.615 multiplication 
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factor for the corrections in oceanic compartments. (Since 
not air but sea water of density 1.027 g/cm^ fills oceanic 
compartments to sea level, only 1.643 g/cm^ additional mass 
density, or 61.5% of 2.670 g/cm^, is needed to complete the 
terrain correction density fill-in.) As one would expect, 
the tables' corrections decrease for a given elevation dif- 
ference with increasing distance from the station, and increase 
for a given zone with increasing elevation differences. Thus 
the elevation difference effectively adjusts the value of the 
numerator of equation (3) , while the zone specification does 
the same to the denominator. 

The effect of terrain upon a gravity value measured on 
the sea floor is more complex than its terrestrial counterpart; 
this is understandable in view of the inclusion of an additional 
medium, the sea water, surrounding the station. This effect 
must be computed in three stages; Figure 13 presents the analy- 
sis in cross-section. 

Had the sea floor (Fig. 13-A) been horizontal, a terrain 
correction would not be required. The existence of a uniform 

O 

slab of 1.027 g/cm J material extending infinitely in all 
directions between station level and sea level can be neglected 
without changing the TC value to be calculated since the topo- 
graphic effect of such a slab of given density and of constant 
thickness equals zero. The first step shown in Figure 13 
simplifies the problem and requires no actual calculations; 
it effectively replaces the sea water in Sector 1 with air and 
decreases the crustal density in Sector 4 to 1.643 g/cm^. The 
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Figure 13. Schematic Representation of the Terrain Correction for a Gravity Station 
on the Sea Floor. (All densities are in g/cm^.) 



resultant mass distribution surrounding the station shown in 
Figure 13-B has an attractional effect (due to terrain) iden- 
tical to that shown in Figure 13-C, i.e., as if the station had 
been located right at the coastline of a layered continent of 
two different densities. 

Step two consists of a "normal" terrain correction as de- 
scribed earlier. Land compartment corrections are sub-totaled 
and added to 61.5% of the sub-total of the corrections for 
oceanic compartments. This step effectively transforms the 
sea water in Sector 2 into rock of density 2.670 g/cm^, and 
removes rock of similar density from Sectors 4 and 5. Sector 
5 is now air (as desired) , but the use of a density of 2.670 
g/cm~* was excessive in Sector 4 , since its density was only 
1.643 g/cm^ after step 1 (Fig. 13-D) . 

Step three corrects this situation by "putting back" 
material of density 1.027 g/cm^, thereby making Sector 4 a 
region of air. This last step is accomplished by subtracting 
from the previously accumulated total correction after step 2 
a correction based on only one elevation (station depth) for 
all landward compartments. It is 38.5% (the density ratio of 
sea water to crustal rock) of the subtotal thus obtained from 
the tables. 

Adding the total terrain corrections to the data refers 
the resulting gravity anomalies to their values had the station 
been located on a flat earth crust (Fig. 13-F) . Table II 
summarizes the steps just described. 
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Sec^o^ ^Designo tion^^^ ^ ^ ^ ^ ^ ^ 
^ ^ ' ' "siatlon 




Step One 


Step Two 


Step Three 


Sector 1 


Decrease 
sea water 
density to 
zero . 


No Change 


No Change. 


Sector 2 


No Change. 


Increase 
water 
density 
to that of 
rock . 


No Change 


Sector 3 


No Change. 


No Change. 


No Change. 


Sector 4 


Decrease rock 
density by 
an amount 
equal to sea 
water density. 


Decrease 
density 
present by 
2.670 g/cm^ . 


Increase 
present densi- 
ty bv 1.027 
g/crtw to yield 
zero effective 
density . 


Sector 5 


No Change 


Decrease 
rock density 
to zero. 


No Change. 



Table II. Summary of Terrain Correction Procedure for a 
Gravity Station on the Sea Floor. (Steps one 
and three are unique to bottom gravimetry.) 
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It would not be possible after step 1 to compute separate 
corrections for Sectors 4 and 5 with density ratio factors of 
0.615 and 1.0, respectively. This is due to the fact that 
the tables were developed with the assumption that the base of 
each land compartment is at station level. 

Most of the maps, charts, and templates required were 
provided by USGS. One such template was constructed by the 
author, however, scaled for use with U. S. Coast and Geodetic 
Survey (USC&GS) Chart No. 5403. The bathymetric information 
on this chart was considerably better than that of others 
available, and was especially helpful for corrections in Zones 
E through I . 

It should be kept in mind that the slope of the bottom in 
Figure 13 is greatly exaggerated. Since bottom slopes in most 
of southern Monterey Bay are on the order of 1:100, many 
stations' cumulative corrections as far out as Zone G totaled 
zero. For this same reason the considerations of steps 1 and 
3 were, in practice, not required. Sector 4 compartment ele- 
vations could never exceed station depth and were always 
negligibly smaller within the closest several zones , with the 
exception of Station 65 (located on the eastern end of the 
thalweg of the Monterey Submarine Canyon). Here, a simple 
multiplication by the 0.615 density ratio factor sufficed for 
the Sector 4 corrections. As this number is usually used to 

transform water in Sector 2 into rock, it similarly transforms 
rock in Sector 4 into water. This simplification was possible 
because Sector 5 corrections for Station 65 are nill (since 
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the elevations of the area surrounding Moss Landing are so 
small) . The fact that the Sierra Nevadas are just beyond 

the 100-mile radius of the most distant zone (Zone 0) and 

\ 

that the San Juaquin Valley is so uniformly flat and close to 
sea level also served to expedite calculations somewhat. 

A realistic assumption regarding possible errors in 
estimating compartment elevations is +0.02 mgal per zone. 

If errors of this magnitude and similar sign existed in each 
of the 15 zones , the total terrain correction error would be 
+0.30 mgal for each station. 

The foregoing analysis treats the terrain correction in- 
dependently of the earlier Bouguer correction. This is valid 
since the Bouguer correction assumes a flat station environ- 
ment (Heiskanen and Vening Meinesz, 1958). Independent analy- 
sis is further justified by considering a change in the order 
of making the corrections: applying a pure terrain correction 

to the free-air anomaly, followed by the Bouguer correction, 
should result in an identical CBA. Regardless of order, the 
free-air correction only looks at vertical distance, the 
Bouguer correction deals with vertical mass attractions, and 
the terrain correction accounts for quasi-horizontal mass 
attractions. 

H. CURVATURE CORRECTION 

The curvature correction (CC) may be thought of as a 
terrain correction of macroscale proportions since it accounts 
for the curvature of the earth. (Recall that the assumed 
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juguer slab was infinitely flat.) For shallow bottom stations, 
lis correction is negative and is given by equation (14) 

3. L. Robbins, personal communication, 1972) : 

CC = 4.462 x 10 -4 H - 3.282 x 10“ 8 H 2 + 1.270 x 10 -15 H 3 

(14) 

lere CC is in milligals and H is height in feet above sea 
;vel . (Oceanic depths are entered as negative numbers.) 

When the terrain and curvature corrections are applied to 
le simple Bouguer anomaly, the result is the complete 
)uguer anomaly (CBA) . Its value may be stated in equation 
jrm: 

CBA = g Q - g t + ET + D - FAC + BC + TC - CC (15) 

:, from equation (12) : 

CBA = SBA + TC - CC . (16) 

le CBA is the difference between observed gravity and the 
ilue expected had the station been located on a flat solid 
irface at sea level, unaffected by astronomical attractions. 
iis final gravity anomaly refines the geological implications 
: the SBA and represents a useful tool that geologists can 
ill upon to infer substructure mass distribution. 

. SHORELINE SURVEY 

Data reduction for the ten stations ashore proceeded in a 
ishion similar to that of the bottom survey, but with five 
cceptions. First, since calibration factors for the con- 
irsion of counter units to milligals are unique to each 
irticular gravimeter, the equation for g Q for the Model G-17B 
:avimeter was not the same as equation (1) . Second, free-air 
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corrections were added (rather than subtracted) , since on 
land, mean sea level in general is closer to the earth's 
center than station elevation. Third, Bouguer corrections 
were subtracted (rather than added) , because part of the 
observed attraction was due to the material between station 
level and sea level. Thus, while underwater Bouguer correc- 
tions are composed of three parts, terrestrial Bouguer correc- 
tions need compensate for only a single attraction. Fourth, 
the three-step terrain correction analysis does not apply. 
Finally, since all elevations were within 10 ft of mean sea 
level, curvature corrections were not required. (A station 
elevation of 22 ft above or below sea level is required to 
produce a curvature correction of 0.01 mgal . ) The equation 
summarizing land station data reduction is then: 

CBA = g Q - g t + ET + D + FAC - BC + TC (17) 

Estimated error values for land station corrections are 
+ 0.02 mgal for g c i +0.07 mgal for g^. (corresponding to 
positioning accurate to within 1/20 n mile) , + 0.06 mgal/ft 
for station elevation, and, again, +0.30 mgal for the terrain 
correction elevation estimations . 

J. DATA PRESENTATION 

The values of the "mass-adjusted" free-air anomaly and 
complete Bouguer anomaly for both shoreline and bottom stations 
were recorded on charts next to appropriate station locations 
and scalar analyses were performed by hand with isolines 
spaced at a 5-mgal interval (Fig. 14 and 15) . Table III 
presents the values of these and two other gravity anomalies 
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of interest for the 92 stations; values of observed and 



theoretical gravity as well as values of the various 
tions are presented in Appendix A. 



correc- 
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Figure 14. "Mass-Adjusted" Free-Air Anomaly Map of Southern 
Monterey Bay. (Contour values are in milligals.) 
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Figure 15. Complete Bouguer Anomaly Map of Southern Monterey 
Bay. (Contour values are in milligals.) 
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TABLF III. DATA PRESENTATION. (VALUES ARE IN MILLIGALS.) 
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VI. DATA ANALYSIS AND GEOLOGICAL INTERPRETATION 



Figure 15 is in general overall agreement with the 
regional trend of CBA values on adjacent land areas as deter- 
mined by Fairborn (1963) , Sieck (1964) , Bishop and Chapman 
(1967) , and Ivey (1969) . The shape of the +10 and +15 mgal 
isolines is in excellent accord with the shape of the granitic 
contact north and east of Pt. Pinos (Fig. 3). Additionally, 
there is a discernible gradient increase from the +10 to -5 
mgal isolines in the southernmost region of the Bay; this 
coincides with the Tularcitos Fault Zone. The previously 
assumed existence of -30 mgal isolines seaward of the Marina 
area (Chapman and Bishop, 1967) was due to insufficient station 
coverage along the coast and no offshore coverage; the -25 
mgal isoline in Figure 15 is confirmed by a CBA value of 
-25.620 at Station H. Monterey Peninsula CBA values were 
found to be slightly higher than those presented on the 1967 
Santa Cruz Sheet, but, again, bottom gravimetry CBA values 
were verified by the independent land survey (Stations A, B, 

C, and D) . Of far greater geological significance than the 
actual CBA value associated with each isoline is the general 
shape and trend of these lines, an aspect of Figure 15 that 
all investigators concerned could easily accept. 

Possible error magnitudes are summarized in Table IV. 

These estimates are considered to be maximum values; it is 
felt that the final CBA values of the author's survey are 
accurate to v/ithin 1 mgal. 
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It is readily apparent from Figure 15 that the north- 
eastern portion of the survey area is characterized by a 
deficiency of mass density relative to the southwestern 
region. This deficiency represents an increase in depth to 
the top of the dense granitic basement from southwest to north- 
east, and is in very good agreement with seismic reflection 
profiles presented by Greene in 1970 (Fig. 16 and 17) . Since 
a single geophysical method seldom assures a unique solution, 
this confirmation of previous seismic reflection interpretation 
is an important result of the author's study. The CBA pro- 
files in the figures lack the detail (particularly in the 
Tularcitos Fault Zone) of the corresponding seismic profiles; 
this is most probably due to the grid size of the station net- 
work. Additionally, there may exist an insufficient density 
contrast, although this possibility is remote in view of the 
apparent displacement and the density difference between the 
Monterey Formation and granodiorite below Monterey Bay. 

Figures 16 and 17 indicate an average complete Bouguer 
gravity anomaly decrease of 5 mgal for each 520 ft increase in 
depth to the granitic basement structure. Since the 11 n mile 
extent of profile A-B and corresponding 3700 ft drop in the 
level of the basement represent a slope at the interface be- 
tween the Monterey Formation and the granodiorite of only 
about 3°, we can approximate the difference in attractive 

forces at two points along that profile whose CBA values differ 
by 5 mgal by the attraction of a horizontally infinite slab 
520 ft thick whose mass density is equal to the density 
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Figure 16. Comparison of CBA and Depth of Granite Substructure 
as Determined by Seismic Reflection for Profile A-B 
(Upper profile after Greene, 1970). The location 
of Profile A-B is shown in Figure 15. 
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Figure 17. Comparison of CBA and Depth of Granite Sub- 
structure as Determined by Seismic Reflection 
for Profile C-D (Upper profile after Greene, 1970). 
The location of Profile C-D is shown in Figure 15. 
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difference between the two rock types. So doing permits 
equating a theoretical attraction of 2lT G p c h (where h is 
520 ft and p c is the density contrast) with the associated 
5 mgal inferred value, and yields a density contrast of 0.75 
g/cm^. Since the slope of the igneous/sedimentary rock inter- 
face is fairly constant along the profile, the values of sub- 
layer terrain irregularity corrections at two such locations 
will be close and their difference can therefore be neglected. 
This analysis assumes, then, that the relation between CBA and 
depth to basement is due solely to the magnitude of the 
density deficiency existing in the upper formation; the 
attraction at the point farther toward the northeast will be 
less since a thicker layer of the (lighter) Monterey Formation 
exists below. 

A fairly accepted value for the density of the granodiorite 
is 2.75 g/crn^ , but the density of the in situ Monterey Forma- 
tion is not accurately known. While this formation consists 
of multiple members, the Miocene marine shale predominates; 
its water-saturated density was determined experimentally by 
the author to be 1.85 g/cm^ (true in situ density is higher 
since the lab samples were not under pressure) . Sieck (1964) 
found the dry density of Monterey shale to vary from 1.41 to 
2.10 g/cm 3 , with an average of 1.80 g/cm^, and his measured 
density of Monterey sandstone averaged 2.10 g/crn-^ . Exami- 
nation of various references to shale densities in general 
indicates that it is a highly variable figure (depending on 
age, location, and many other factors); values for diatomaceous 
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shales in the San Jauquin Valley are between 0.9 and 1.1 g/cm 3 
while some water saturated shale types can be as dense as 
3.21 g/cm 3 (Jakosky, 1950). In view of these many factors, 
the 2.00 g/cm 3 implied density of the Monterey Formation (2.75 
g/cm 3 granodiorite minus the gravity survey's calculated 0.75 
g/cm 3 density contrast) would appear quite acceptable as a 
true value. 

Having arrived at densities of 2.00 and 2.75 g/cm 3 for 
the two major geological units in the area, it is now appropri 
ate to consider the validity of using 2.67 g/cm 3 for the 
assumed average crustal Bouguer density. Since the Monterey 
Formation is less than 1 km thick, while the granodiori te sub- 
structure is probably several times thicker, the average 
crustal density in the vicinity of southern Monterey Bay is 
much closer to 2.75 than to 2.00 g/cm 3 . Thus, while the 
frequently used average crustal density of 2.67 g/cm 3 can be 
unrealistic in many regions of the world, it is probably not 
the case here. 

Throughout the course of this research it was anticipated 
that isoline positions over the deep Monterey Submarine Canyon 
could be inferred from the results of the bottom surveys in 1 
the northern and southern halves of Monterey Bay. Both sets 
of data were obtained with the same equipment in an identical 
manner and reduced in a comparable fashion, so there is no 
reason to doubt the validity of making such an inference. 
Figure 18 is the result of combining the data from Cronyn ' s 
(1973) survey with that of the author. 
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Figure 18. Composite Gravity Map of Monterey Bay (Northern 
Monterey Bay contours after Cronyn, 19 73) . 
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The fact that the southernmost positions of Cronyn's known 
0,-5, -10, and -15 mgal isolines must bend to meet the northern- 
most positions of those of the author supports the theory that 
a fault exists along the axis of the Monterey Submarine Canyon. 
The difference in lithologies on opposite sides of the Canyon 
led Martin and Emery (1967) to recognize the so-called Monterey 
Fault parallel with the Canyon and offset at numerous locations 
by transverse faults of different magnitudes. Although the 
fault's suggested orientation is at an angle to the trend of 
almost all of the many fault zones in central California, the 
proposed situation is seen elsewhere, as in the case of the 
Santa Ynez Fault region north of Los Angeles, where the strike 
of the San Andreas Fault changes. Further evidence of the 
existence of the Monterey Fault based on seismic reflection 
tracks across the Canyon has been presented by Greene (19 70) . 

Although some of the advocates of the Monterey Fault 
suggest that the motion along the fault is dip-slip with the 
northern side having descended relative to the southern side, 
this is not uniquely indicated by the gravity data. On both 
sides of the Canyon (in the eastern portion of the Bay) the 
CBA isolines run north-south with the higher values seaward 
(Fig. 18 and 19-A) . If the motion of the Monterey Fault is, 
in fact, primarily dip-slip, the shape of the CBA field over 
the Canyon would suggest a relative motion such that the 
southern side has descended rather than the northern side, as 
shown in Figure 19-B. This is supported by the fact that the 
motion of the northeastern side of the Tularcitos Fault is 
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Figure 19. Two Possible Structural Explanations for the 
Inferred CBA Field in the Vicinity of the 
Monterey Submarine Canyon. 
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downward (Greene, 1970). If, on the other hand, the fault 
motion is primarily right-lateral strike-slip, the inferred 
CBA isolines over the Canyon may be explained as in Figure 
19-C, wherein the northern side has moved east relative to the 
southern side. Specifying which of the two types of relative 
motion has occurred is not possible from the gravity infor- 
mation alone; quite possibly a combination of both types of 
motion has occurred, causing an oblique fault. 

It may well be that the intersection of the Monterey Fault 
and the Palo Colorado - San Gregorio system farther to the 
west corresponds to the area where a concentration of recent 
seismic epicenters is located. More geophysical and ocean- 
ographic study of this region is needed before an accurate 
structural model can be developed. 
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VII. FUTURE WORK 



It is recommended that CA-259 (the Monterey County Airport 
station) of the state gravity base station network be re- 
occupied in accordance with procedures outlined by Chapman 
(1966) to re-establish a correct value. Local construction 
there will be completed soon. 

It is suggested that additional gravity measurements be 
made in the vicinity of the Monterey Submarine Canyon. Where 
a large portion of the coverage needed might be precluded 
(with a bottom gravimeter) by excessive depths and/or slopes , 
surface ship gravimetry could be called upon to complete this 
task. One such survey consisting of several canyon crossings 
was conducted by USGS in November, 1972, aboard the USNS 
BARTLETT, but at this writing the results are not yet available 
for inclusion here. They will eventually help greatly in re- 
moving the question marks in Figure 18. 

At this writing the gravimeter and auxiliary equipment 
used for the Monterey Bay surveys is back aboard the R/V 
ACANIA. The author has thoroughly familiarized two new NPS 
students with the entire operation, and surveys are underway 
to extend shallow water coverage in two areas, one from Pt. 

Lobos to Pt. Sur and the other from Santa Cruz northwest to 
beyond Davenport. The goal of the project sponsors is to tie 
in these surveys with those conducted by the author and co- 
workers, so that a continuous picture of the off-shore Bouguer 
anomaly will become known. Hopefully this goal will be achieved. 
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APPENDIX A (Continued) 



SHORELINE SURVEY 
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3.81 
3.78 
3.21 
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2.84 
2.50 
2.12 
2.09 
1.97 
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COMPUTER PROGRAMS 



CCMPUTER PROGRAM FOR SHORELINE STATIONS 



IMPLICIT REAL; 3 (A-H.O-Z) 

DIMENSION THEO(IO) , E L E V ( 10) , FAC ( 10 ) , F AA ( 1 0 ) » BC ( 10 ) , SB A 
1( 10) , C B A ( 1 0 ) , T L A ( 1 0 ) , A RGA ( 10) . ARG5 ( 10 ) , DE GR E E ( 1 0 ) , GO 8 ( 
1 10 ) , ET ( 10) ,TC( 10 ) 

READ (5 *500) (GC3( I ) .1 = 1. 10) 

READ (5.1000) (DEGREE ( I ) , I =1 ,10) 

READ (5,2000) ( EL E V ( I ) , I = 1 , 1 0 ) 

READ (5,3000) ( E T ( I ) , I = 1 , 1 0 ) 

READ ( 5, 3500) (TC( I) ,1=1,10) 

500 FORMAT (8 F10.3) 

1000 FORMAT ( 8F l 0.7) 

2000 FORMAT ( 1 0 C 3 . 1 ) 

3000 FORMAT (10F4. 2 ) 

35C0 FORMAT ( 10F4. 2) 

DO 4000 1= 1. 10 

TLA ( I ) =3.14159 DO- ( DEGREE ( 1)7180.0 DO) 

ARGA ( I ) =C SI N( T LA ( I ) ) 

ARGB ( I ) = DS ! N { 2 . 0 DO ' TL A ( I ) ) 

THEO ( I ) =978049.000;' (1 . 0D0 + ( 0 .00 52384 D0 : ' ( (ARGAU ) )**2) ) 
1-C. 00000 59D0" ( ( A R G 3 ( I ) ) 2 ) ) 

FAC ( I ) =0 .09405 DO: ( EL EV ( I ) ) 

F A A ( I) =GG8 1 I ) -THE 0 ( I ) +ET ( I ) +FAC ( I ) 

BC ( I ) = ( 2 . 0 DO'- - 3 . 1 4 1 5 9D O' • 6 . 6 7D 0 30.48D0'' 2. 67 DO"- ( ELE V( I ) ) 
1 ) /100000 .Or'O 
SR A ( I ) = F A A ( I ) - BC (1 ) 

4000 CBA( I )=S3A( I ) +TC ( I ) 

DO 5000 1=1,10 

WRITE ( 6,6000) I ,GQ8 ( I ) , THE 0 ( I ) , ET ( I ) , E L E V ( I ) , FAC ( I ) , F 
1 AA ( I ) , BC ( I ) , S B A ( I ) , TC ( I ).C6A( I ) 

5000 CONTINUE 

6 000 FORMAT ( / // , 20X , I 2 , 3 X , F 1 0 . 3 , 3X , F 1 0 . 3 , 3X , F4 . 2 , 3X , F3 . 1 , 3X 
1,F6.3,3X,F7.3,3X,F6.3,3X,F7.3,3X,F4.2,3X,F7.3) 

STOP 

END 
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COMPUTER PROGRAM FOR UNDERWATER STATIONS 



500 

501 
1000 
1200 
1201 
1202 
1203 
2500 
4444 

7000 

7001 

7004 

7005 

7006 



IMPLICIT REAL I 8 (A-H.O-Z) 

DIMENS ION DATE (90 ) . GMT (90 ) , L AT DEG ( 90 ) , L ON DEG ( 90 ) , L A TMI 
IN (90) . L A TSEC (90) , LON MI N(90) , LONS EC (90) , COUNT R (90 ) • GR EL 
1(90), V A L GR V ( 90 ) .TERCOR (90) , BUGA ( 90) ,GVTI DE (90) , BUGAET ( 
190), DECK (90 ) .BOTT (90) , D I F F ( 90 ) , DEPTH ( 90 ) , T IDE HT ( 90 ) , FA 
1 C ( 90) ,FAA( 90) , BC (90) , S3A(90) , DEGREE (90) , TLA (90) ,ARGA(9 
10), ARGB ( 90 ) , THEO( 90) ,CURV( 90) ,CBA ( 90) , TFAC (90) ,TF AA( 90 
1 ) 



READ (5,500) 
READ (5,501) 
READ (5,1000) 
READ ( 5. 1200) 
READ (5,1201) 
READ ( 5,1202) 
READ (5, 1203) 
READ (5,2500) 
READ ( 5,4444) 
READ (5.7000) 
RFAD (5.7001) 
READ ( 5, 7004) 
READ (5,7005) 
RFAD (5,7006) 
FORMAT ( 20A4) 
FORMAT (20A4 ) 
FORMA T ( 1 1F7.2 
FORMAT ( 40 A 2 ) 
FORMAT (40A2) 



( DATE ( I ) ,1=1,82) 

( GMT ( I ), 1 = 1,82) 

( COJNTR( I ), 1 = 1 ,82 ) 
( L ATM INC I ) , I =1 , 32 ) 
( L ATS EC ( I ), 1 = 1,32) 
( LONMIN ( I ) ,1=1,82) 
( LON S EC ( I ) , I = 1 , 82) 
( D E G R E E ( I) ,1 = 1,82) 
(GVTIDE(I) ,1=1,82) 
(DECK ( I ) , 1 = 1, 82) 

( BOTT ( I ) , 1=1,82) 

( T I D F HT ( I ) ,1=1,82) 
( T ERCOR ( I ) , 1 = 1,82) 
(CURV(I) ,1=1,82) 



, 3 X) 



FORM AT ( 40A 2 ) 
FORMAT (40A2) 



FORMAT ( 8 FI 0.7) 



FORMAT ( 2 0 F 4 . 2 ) 

FORMAT (20F4.1 ) 

FORMAT ( 16F5. 1 ) 

FORMAT ( 20 F4. 1 ) 

FORMAT (20F4. 2) 

FORMAT ( 16F5. 3) 

DO 7500 1= 1,82 

LATDEG ( I ) =36 
LONDEG( I ) = 121 

TL A ( I ) =3 .1415900 ( DEGREE ( I )/ 18 0.0 DO) 

ARGA ( I ) =DSI i\(TLA( I ) ) 

ARGB ( I ) = DSI N ( 2 . OD O' TLA ( I ) ) 

TF'EO ( I ) =973049 .ODO- ( 1 . 0 D0+ ( 0 . 00 52884D0 » ( ( A RG A ( I ) )**2) ) 
1-0. CG0005900 ( ( A RG B ( I ) )** 2 ) ) 

VALGRVl I ) = 9793 9 1 . 7000 + ( ( COUN TR ( I ) -33 23. 0 500 0 ) '-> 1 . 039850 

IDO) 

B UGA ( I ) = VA LG RV ( I ) - TH EO ( I ) 



BUGAET ( I ) = BUGA ( I ) +GVT IDE ( I ) 

D I F F ( I 1 = BOTT ( I ) — DE CK ( I ) 

DEPTH ( I ) = 0. 7700 0" 01 F c ( I ) 

F AC ( I ) =-0. 0940 60D0 M DEPTH ( I ) -T I DEHT ( I ) ) 

FAA( I ) = 3UG AET( I ) +FAC( I ) 

TFAC ( I ) =-0.094 060 0 (DE°TH( I )- T IDEHT(I ) ) + 2 . ODO* 3. 1 4159D 
1 0-: 6 ,67D0v- 1 .027 D0 : '30 ,48D0~-DEPTH ( I ) / 1 00000 . OD 0+ 2. ODO'^S. 1 
141590 0* 6.670 0*1 . 02 7D0? 3 0. 4 8D0* ( D E PTH ( I ) -T I DE HT ( I ) ) / 1 00 
1000 .ODO 

T F A A ( I ) = BUGAET ( I ) +T FAC ( I ) 

BC( 1 ) = ( ( 2. ODO* 3. 141 59 DO - ' 6. 6700? 1 . 027D0 •' 30 . 48 DO" DE PTH ( I 
1) )/ 100000. ODO) +( (2.0D0- ; -3.14159D0--6.67D0*2.67D0 ? -30.48C0 
IT ( DEPTH ( I )-TIDEHT( I ) ) )/l 00000. ODO ) 

SB A ( I ) = FAA( I ) + RC ( I ) 

I F ( GVT I D E ( I ) .59.0.0 ) GVT I DE ( 13 = 1. E~ 10 
7 500 C 3 A ( I ) =SBA( I )+ TERCOR (I ) + CURV(I) 

STOP 

END 



(OUTPUT FORMAT STATEMENTS HAVE BEEN OMITTED.) 
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